MONIC NON-COMMUTATIVE ORTHOGONAL POLYNOMIALS 



MICHAEL ANSHELEVICH 

Abstract. Among all states on the algebra of non-commutative polynomials, we characterize the 
ones that have monic orthogonal polynomials. The characterizations involve recursion relations, 
Hankel-type determinants, and a representation as a joint distribution of operators on a Fock space. 



1. Introduction 

For a measure on M all of whose moments are finite, there are two standard ways to normalize 
the polynomials orthogonal with respect to fi. One can take the polynomials to be monic, 

Pn{x.) = + lower order terms. 

Or one can take them to be orthonormaJ, 

/ Pn{xY d^{x) = 1. 

For a measure ji on R'^, the situation is more subtle. One can always orthogonalize the subspaces of 
polynomials of different total degree (so that one gets a family of pseudo-orthogonal polynomials). 
The most common approach is to work directly with these subspaces, without producing individual 
orthogonal polynomials; see, for example [DXOl]. One can also further orthogonalize the polyno- 
mials of the same total degree, for example to make them orthonormal [BCJ05, BC04]; however 
this requires a choice of an order on the monomials of the same degree, and there is no canonical 
choice of such order. The third approach, and (it is easy to see) the only one that will produce monic 
orthogonal polynomials, is to require that the pseudo-orthogonal polynomials obtained in the first 
step already be orthogonal. The price one pays is that this can be done only for some measures /i. 

This paper is about orthogonal polynomials in non-commuting variables. There is no difficulty with 
the definition. The usual orthogonal polynomials are obtained by starting with a measure fi on W^, 
thinking of X2, . . • , x^] as a vector space with the (pre-)inner product 

(P,Q)= / P(x)g(x)d/i(x), 

and applying the Gram-Schmidt procedure to the monomials {^Xu{i)Xu{2) ■ ■ -Xuin)}- In the non- 
commutative case, one starts directly with a positive linear functional (state) ip on the algebra of non- 
commutative polynomials M(xi, X2, ■ ■ ■ , Xn), and orthogonalizes the monomials in non-commuting 
variables with respect to the inner product 

(P,Q) = (^[P*(x)Q(x)]. 
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In particular, in [Ans06], I showed that monic polynomials in non-commuting variables orthogo- 
nal with respect to a faithful state satisfy a recursion relation, as orthogonal polynomials should. 
The question remained: which states have monic orthogonal polynomials? It is the first question 
answered in Theorem 2 of this paper, where moment conditions and Fock space representations of 
such states are provided. The second question, asked by the referee, was whether one needs the 
faithfulness condition. In this paper, that condition is removed using a new technique, namely a 
representation of the state as a joint distribution of some operators. Such a representation is closely 
related to a combinatorial way of representing moments as sums over lattice paths [Ra80], but in 
the multivariate case I find the operator formulation more useful. 

An important class of states that have monic orthogonal polynomials are the free Meixner states, 
whose study I initiated in [Ans06] . Applications of the techniques from the current paper to free 
Meixner states are considered in a companion paper [Ans07]. 

2. Preliminaries 

2.1. Polynomials. Let x = (xi, 0:2, . . . , x^) be a d-tuple of non-commuting variables. Let 

R(x) = R{xi,X2, ...,Xd) 
be all the polynomials with real coefficients in these variables. Multi-indices are elements 

■u G {1, . . . , d}'' 

for A; > 0; for |m| = denote u by 0. Monomials in non-commuting variables (xi, . . . , Xd) are 
indexed by such multi-indices: 

Note that our use of the term "multi-index" is different from the usual one, which is more suited for 
indexing monomials in commuting variables. 

For two multi-indices u, v, denote by {u, v) their concatenation. For u with \u\ — k, denote 
Define an involution on M(x) via the R-linear extension of 

A monic polynomial family in x is a family {P3(x)} indexed by all multi-indices 

00 

[j{ue{l,...,d}'} 

k=l 

(with Pi/, — 1 being understood) such that 

Pjj(x) — Xii + lower-order terms. 
Note that Pi ^ P{ii)op in general. 
Definition 1. A state on M(x) is a functional 

(/9 : R(a;i,a;2, . . . ,Xd) 
that is linear, compatible with the ^-operation, that is for any P, 

cp[P]^^ [P*] , 
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unital, that is <^ [1] = 1, and positive, that is for any P, 

[P*P] > 0. 

A state is faithful if in the preceding equation, the equality holds only for P = 0. Unless noted 
otherwise, the states in this paper are not assumed to be faithful. 

The numbers (p [xn] are called the moments of (p. 
Remark 1. A state (p induces the pre-inner product 

{P,Q)^^<p[P*Q]^{Q,P)^ 

and the seminorm 

Throughout the paper, we will typically drop </? from the notation, and denote the inner product and 
norm it induces simply by (■,■), ||-||. 

Even if a state is not faithful, it is easy to check that the Cauchy-Schwartz inequality 

\{P,Q)\<\\P\\ WQW 

still holds. In particular, if ||P|| = 0, then for any Q, (P, Q) — 0, and if ||P — P'|| =0, then for any 

g, {p,Q) = {p',Q). 

Definition 2. A state has a manic orthogonal polynomial system, or MOPS, if for any multi-index 
u, there is a monic polynomial P^ with leading term x^j, such that these polynomials are orthogonal 
with respect to </?, that is, 

(Pff, Pp) = 

for u^v. 

Note that the same abbreviation is used in [DES] to denote a class of multivariate orthogonal poly- 
nomials systems, which is different from ours. 

Lemma 1. Let {Qu} cind {P^} he two monic polynomial families, and suppose that {Qa} are 
orthogonal with respect to a state </?. The following are equivalent: 

(a) {Pa} are pseudo-orthogonal with respect to (p, which means that 

{Pa, Pp) = 

whenever \u\ ^ \v\. 

(b) For each u, Qa and Pa are equal in L'^{(p), that is, 

\\Qa-Pu\\^0. 

If either of these statements holds, then in fact {Pa} are orthogonal. 

Proof (a) => (b). Since Qa and Pa are monic, Qa — Pa has degree less than \u\ . Since both families 
are pseudo-orthogonal, Qa — Pa is orthogonal to all polynomials of degree less than \u\. Therefore 
it has norm zero. 

(b) =^ (a). Using Remark l,forv^u, {Pa, Pa) — {Qa, Pa) — {Qa, Qa) — 0. Note that this in fact 
implies that {Pa} are orthogonal. □ 

Lemma 2. Let ip be a state with MOPS {Pa}- Suppose that for some u, \\Pa\\ = 0. Then for any i, 

II ^M) II =0. 
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Proof. Since {Pi;} are monic, they form a basis for the vector space of polynomials. So we may 
write 

|^;|<|«| 

Then 

|i;|<|3| 

= (a;i-PM)(x),Pj:(x)) - ^ a^; (P(j,3)(x), P^;(x)) = 

|ff|<|«| 

because of orthogonality and Remark 1. □ 

3. MONIC ORTHOGONAL POLYNOMIALS STATES 

The goal of this paper is to prove Theorem 2, which provides a number of equivalent conditions 
describing a class of states. Since these conditions come from quite different frameworks, we first 
describe two constructions. Among other things, these results now apply to not-necessarily-faithful 
states, thus answering a question of the referee of [Ans06]. 

Notation 3. Put an order on all multi-indices that is compatible with degree and otherwise arbitrary, 
for example the lexicographic one. Let \u\ = n. Construct dil + d + cP + . . . + d"^^ + 1 = + 1- 
dimensional square matrix A^, whose rows and columns are labeled by all the multi-indices of 
length less than n and u, and whose {v, w) entry is {xp, x^) . Denote /i^j = det A^x, and f)„ the 
determinant of the matrix A^t with the w'th row and column removed (so that it depends on n but 
not on u). These are multivariate versions of Hankel determinants, (compare with [BC04]). Note 
that hn, i)n do not depend on the chosen order on multi-indices. Let A^u(x) be the matrix with 
the u'th row changed, so that the {u, w) entry is x^. Finally, let be the determinant of A^^ with 
the u'th row changed, so that the {u, w) entry is {x^, x^). In particular, h^t^Q — h^. 

3.1. Fock space construction. Let H — C^, with the canonical orthonormal basis ei, 62, . . . , e^. 
Define the (algebraic) full Fock space of H to be 

A;=0 

Equivalently, J^aig(^) is the vector space of non-commutative polynomials in ei, 62, . . . , e^. Fol- 
lowing convention, we will denote the generating vector in 7Y®° = C by Q instead of 1. 

For each A; > 1 let C^*^) be an operator 

We identify each C^*^) with its x d^ matrix in the basis {e„(i) ® . . . eu(fc) }• Assume that for each 
k, C^^^ is diagonal and C^^'^ > 0. 

For i — 1,2, . . . ,d, define af and a~ to be the usual (left) free creation and annihilation operators, 

af (en{i) ® e„(2) ® . . . ® e„(fc)) = Cj e„(i) (g) e„(2) ® . . . ® e„(fc), 

a," (eu(i) (g) e„(2) ® . . . ® eu{k}) = (e^, eu{i)) eu{2) ® • • • ® eu{k)- 
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For each i — 1,2, . . . ,d and each A; > 0, let 7^^ ' be an operator 



Assume that each satisfies 



(1) {r^'^y ICc ^ }CcT^'\ 

where /Cc is the operator in equation (2) below, and A* is the transpose. We will denote by % and C 
the operators on J^aigCH) acting as T^''^ and C'^'^^ on each component. Finally, denote 

Note that a~^l = a~^l = 0. 

On each H®'' one has the usual inner product (•, •) induced from H. Define a new pre-inner product 
(■, by using the non-negative kernel 

(2) /Cc = (/ ® / . . . ® / ® C(^)) . . . (/ ® / ® C('^-2))(/ ® C^k-^^)C^k'^ 
where / denotes the identity matrix of the appropriate size. That is. 

Put this pre-inner product on each component of J^a\g{H) and define the components to be orthog- 
onal among themselves. Factor out the vectors in J^^^giTi) of norm zero, and denote the completion 
of the factor with respect to this inner product TciV). 

Lemma 3. Let S — {rj & J-3Xg{Ti.) \ W^lWc ~ 0}, so that J-e{Ti) is the completion o/.Faig(7i) / 5". 

(a) S = Span (e„(i) ® ... e^u[k)\k > 0, ||e„(i) . . . e„(fc)||^ = O) = ker/Cc. 

(b) The operators aj' , %, factor through to Tc {'H). 

Proof, (a). Let rj E S, 

n 

fc=0 \u\=k 

Since each C^^'> is diagonal, we can denote C(e„(i) ® . . . e„(fe)) = C^e„(i) ® . . . e^^k) and 

^c(eu(i) ® . . . e„(fc)) = /Cc,^e„(i) ® . . . e„(fc). 

Then 

In n 

= WvWc ^ ( X] XI <8) ■ ■ ■ e„(fe), /Cc ® ■ ■ ■ e^(j) 

\ fe=0 |u|=fe j=0 \v\=j 

n n \ 

?7we„(i) (g) . . . eu(A;), ^ ^ ICcvVve^ii) ® • • • / = X] ifu^c,u- 

, fc=0 |3|=fe j=0 |v|=i / u 

Since /Cc,« > 0, it follows that 

ry- ^ =^ /Cc,n = ^ ||e„(i) (g) . . . e„(fc) ||^ = 

and also that K-cifl) = 0- 
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(b). By part (a), it suffices to consider the actions on e„(i) (g) . . . eu{k) € S. Then 

Kicat{^u{\) ® ■■■ eu{k)) = A:c)C(ei <8) e„(i) . . . e„(jfc)) 
= C(i,it){I <8) A:c)(ei (g) e„(i) O . . . e„(fc)) = C(j,5)a/"A:c(e„(i) (g) . . . e„(fe)) = 0. 
Also /Cc^ = ^*^c by definition of %, so that 7^(ker /Cc) C ker /Cc. Finally, it is easy to check that 

(3) }Ccal = lCca~C = a^{I ® lCc)C = a^^c 

and so a~(ker /Cc) C ker /Cc. □ 

For each i, define an operator A'j on Tc {H) with dense domain J^^igiTi-) /Shy 

Xi = al + % + a^. 
Proposition 1. Each Xi is a symmetric operator. 
Proof. Using equation (3), for any u, v, 

{at (e„(i) (8) ... (8) e„(fe)) , e^(o) (8) . . . O e^(jfc)>^ 
= (cj (8) e„(i) (8) ... (8) e„(jk), e„(o) (8 ... (8) e„(fe)>^ 
= (ej (8) e„(i) (8) ... (8) e„(jk), /Cc (e^(o) (8 ... (8 e^(jt))) 
= (e^,(i) (8) ... (8) eu{k), (h^c (ei,(o) 8) ... (8) e^(fc))) 
= (e„{i) (8) . . . <8) e„(fc), /Cca," (e^(o) (8) . . . <8) e^,(fc))) 
= (e„(i) 8) ... (8) e„(fc), (et,,(o) (8) ... (8) (^v{k)))c > 
so with respect to the C-inner product, + is symmetric. Similarly, 
{% (en(i) 8) ... 8) e„(fc)) , e^(i) (8) ... (8) e^,(fc))^ 

= (e„{i) (8) ... 8) e„(fc)) , /Cc (e^,(i) (8) ... (8) e^(fc)) ) 
= (e«(i) ... 8) eu(k),T/}Cc (e„(i) (8) ... e„(fc)) ) 
= (e«(i) <8) . . . <8) eu{k),}Cc% (e„(i) . . . <8) e„(fe))> 
= (e„(i) (8 . . . <8) e„(fe),'7; (e^(i) (8) ... e^(jk))>^ , 
so % is symmetric. □ 

Definition 4. For any choice of the matrices C^''^ and as above, the corresponding (Fock) state 
if — <pc,{ri} on ]R(x) is defined by 

(fi [P{xi, X2,..., Xd)] = P{Xi, X2,..., Xd)^)(, . 
Lemma 4. </? satisfies all the properties in Definition 1. 

Proof. Since </? is a vector state corresponding to the vector Q and = 1, <^ is linear, unital, and 
positive. Finally, since each Xi is symmetric, 

if \Xu(l) ■ ■ ■ Xu(k)\ = '^m(I) • • • 'V«(fc)^)c ^ (-^"(fc) • ■ • '^u(l)^, 

= (fi, Xu(k) ■ ■ ■ Xu(i)^)^ = If [{xu{i) ■ ■ ■ Xu{k))*'\ ■ n 

Lemma 5. For any C, {%} as above, the state fc,{ri} has a MOPS {Pit}, which satisfy 

(4) Pu{Xi, X2,..., Xd)fl. = e„(i) 8) ... 8) eu{\a\). 
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Proof. For each u, 

X^iQ, = 6^(1) ... e„(|5|) + r]u 

for some r]u G @fl~^ n®\ So we can recursively construct (non-unique) monic polynomials 
P^(x) which satisfy equation (4). Note that since all Xi are symmetric, -P|(X) really is the adjoint 
of P^(X). Then 

= (P3(X)0, P^;(X)0)g = (e„(i) (g) . . . e„(|3|), e„(i) ... e„(|j7|))^ = 

for u^v. Thus {Pj:} are monic orthogonal polynomials. □ 

In the following theorem, the most interesting characterizations are part (g), which uses only mo- 
ments of </?, and part (e), which provides a way to construct any (/? from matricial data. 

Theorem 2. Let ip be a state on ]R(x). The following are equivalent: 

(a) The state (p has a monic orthogonal polynomial system. 

(b) The polynomials {Py} defined recursively by the Gram-Schmidt relation 

(5) Pu -^u ^ ^ 1 1 7^ 1 1 2 

v: \v\<\u\, IIPsIIt^O II^^II 

are a monic orthogonal polynomial system for ip. 

(c) There is a family of polynomials {Pu} such that if [P^] = 0/or allu ^ % and they satisfy a 
recursion relation 

Xi — Pi + Bi^iif^ti,, 

d 

(6) -^iPu P{i,u) ~l~ ^ ^ Pi,w,uPw ~l~ ^i,u^m 

w=l 

^iPu — P{i,u) '\~ Bi,iS,uPiS '\' ^i,u{l)CiiP(u{2),u{Z),...,u{k)): 

|to| = |it| 

with > and, denoting Sj — {s{j), . . . , s{k)), 

k k 
Bi,s,u Cg. — Bi^^ g Ci^. . 

J=l J=l 

(d) The state ip has a MOPS {Pa} and for any u ^ w, \u\ = |w| = n, 

,q\ I \ \ ^ {^u: Pv) {Pvy^w) 

\') \Xu-,XiS) — /p^ p\ ' 

\v\<n,\\P^Ui) ^ 

SO that the even non-symmetric moments ofip are determined by the rest of its moments. 

(e) The state ip has a Fock space representation (pc,{Ti} as in Definition 4. 

Ifip is faithful, the following two conditions are equivalent to the preceding ones: 

(f) The polynomials | ^ det Mjj(x) | (from Notation 3) are orthogonal with respect to ip. 
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(g) For any u ^ w, \u\ — \w\ — n, 
(8) (X5,x^)= 2^ — — — . 

Proof, (a) =^ (b). Suppose that 93 has a MOPS {Qs}. Proceed by induction on n = \u\. = = 
1. Suppose that (P^^, = for all \w\ , |w| < n. From relation (5), it follows that (P^^, P^) = 
for \v\ < \u\. Using one direction of Lemma 1, we conclude that \\Pg — Q^\\ = 0. Using the 
other direction of Lemma 1, it follows that for any u ^ v with |m| , I^T] < n, {Pjx, Pjj) = 0, and the 
induction hypothesis is satisfied. Therefore {Pn} are a MOPS. 

(b) =^ (d). By Lemma 1, all MOPS for (p have the same inner products, so we might as well use the 
MOPS from part (b) which satisfy the Gram-Schmidt recursion (5). Then for u w, \u\ = \w\ ~ n, 

I \ /p P\_L {•^W) P'ti) {Puj Pv) V~*' (•^u; Pv) (P-vi Pw) 

\Xui Xy}) — X-MT, -fiB/ + / ^ ,, ,,2 T" / ^ 



v: \v\<\w\, WP^UO ll^'^ll v: \v\<^\\P^\\^0 H^'^H 



(9) 



Pv) (Xyj^ Pg 

Xu, Pv) {P-Vy -^li 



IIPII^ IIP 11^ X-MT, -Ts/ 
v,s: \v\,\s\<\w\, \\P^\\^0^\\Ps\\ W^vW limi 



^ (P- P-) 

\v\<n,\\P^\\^0 \ ^' 

(d) =^ (a). Obvious. 

(a) =^ (e). Let 9? be a state with MOPS {Pu}. For each k, define a diagonal matrix C'^*^^ recursively 
via 

(10) [-^sW-^^W] = (^"(1) ® ^c{eu{2) ® . . . ® e„(fe)),C(*^)(e„(i) ® . . . ® e„(fc))> 

= ||e„(i) ® . . . (8) e„(fc)||^ . 

This may not be well-defined if /Cc(e„(2) ® . . . e„(fc)) = 0. However, in that case \\P(u(2),...u(k)) \\ = 
0, which by Lemma 2 implies that if [Pi(x)P^(x)] = 0, so one can take the corresponding entry of 

C^'^^ to be anything, for example zero. Similarly, for each k and for \u\ — \v\ — k, define T^''^ via 
^^^^ 93 [Pi(x)a;iP^r(x)] = (/Cc(e„(i) . . . (g) e„(fc)), 7^(e^,(i) (g) . . . (g) ei,(fe))) 

= (eu(i) ...0 e„(fe), 7;(e^(i) (g) . . . (g) e^,(fc)))^ . 

Again, this may not be well-defined if ||e„(i) (g) . . . (g) e„(A;) ||^ = or equivalently if ||Ps(x) || = 0. 
But in that case cp [Pi(x)a;iP^(x)] = (P^, x^P^) = 0, and one can take the corresponding entry of 
T^'^'' to be anything, for example zero. 

By construction, each C*^^^ is diagonal and non-negative. Also, 

V [Pi(x)x,P^(x)] = [(a;iP„-(x))*P^(x)] = ^ [P;(x)a:,P„-(x)] 
so {}Cc%y = VJCc = lCc%. 

Let {Qa} be a MOPS for </7c,{ri} constructed in Lemma 5. Then for any |{t| — \v\, the equations (11) 
and (10) translate to 
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and 



(13) 



We may assume that both {P^} and {Qu} satisfy, for their respective states, the Gram-Schmidt 
recursions (5). Now we prove, by induction on the degree, that one can take = -Pa, which 
impUes the equality of the states — <pc,{Ti}- Suppose that — for \u\ <n. Then it follows 
from equations (12) and (13) that the two states coincide on polynomials of degree at most 2n + 1. 
But in that case, the Gram-Schmidt recursions for Pq and Qn, for\u\ — n + 1, are identical. 



(e) ^ (c). Let </? = ^c,{%}- Let {-Bi,ti;,u, C^} be the matrix elements of , so that 



C(e„(i) ® ...® e„(jfc)) = Cae„(i) ® ...® eu{k)- 

Note that the conditions on the coefficients in part (c) correspond exactly to the conditions on op- 
erators in Section 3.1. Define the polynomials {Pa} using the recursion (6). We show by induction 
that relation (4) holds. Indeed, 



It follows that for all Uy^$,(p [PsCx)] = 0. 

(c) =^ (a). For this direction only, the arguments of Proposition 3 of [Ans06] apply, and it fol- 
lows that {Pa} are orthogonal. Note that in that proposition (p was assumed to have zero means 
and identity covariance; however, it is easy to modify its proof to get a state even without those 
assumptions. 

Now assume that a state (p is faithful. Then for every n, 7^ 0. Indeed, if some f)„ = 0, then for 

the corresponding matrix, a linear combination of some of its rows is zero. This is equivalent to 
saying that for some polynomial P of degree less than n, {P, x^) = for all 1^1 < n. But in that 
case, ||P|| = 0, which contradicts the faithfulness assumption. 

f)|3| is the leading coefficient of detMff(x), and since it is non-zero, l^detMsj are a monic 

polynomial family. Moreover, for \v\ < \u\, (x^, det M,g) equals L^^g, the determinant of the matrix 
A,g with the w'th row replaced so that the w) entry is (x^, x,g) . Such a row is identical to the iTth 
row, and so the determinant is zero. Thus these polynomials are pseudo-orthogonal. 

(a) (f). One implication follows from Lemma 1. The other one is trivial. 




and 




= -^i(e„(i) ® ...® eu{k)) - %{eu{i) ■ ■ ■ <8) e„(jt)) - C(e„(i) (g) . . . (g) eu{k)) 
= at(e„(i) (8) ... e„(jk)) =6^(8) e„(i) ... e„(fe). 
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(f) 4^ (g). Given Lp, let {Pa} be a pseudo-orthogonal family defined via the Gram-Schmidt recur- 
sion (5). From the calculation (9), it follows that 

^ \ _ / p p \ I \^ {Xu.Py) {Pv,X^) 

\v\<n,\\P^\\^0 ^ ^' 

Thus for pseudo-orthogonal polynomials, condition (7) is equivalent to orthogonality. It remains to 
note that 

{xu, det Mil) = hu,v 

and in particular 

(x^;, det Mj;) = 

so that 

hv_ 

h\v\ 

and condition (8) is exactly the condition (7) for the polynomials | ^ det Mg| . □ 
Example 2. For degree one, parts (d, g) of the theorem say that for j, 

so that the variables are uncorrelated. For degree two, they say that for (i, j) ^ (t, s), and assuming 
for simplicity that (p [xi] — for all i. 



7^ det M, 



ifi [xiXjXgXt] = (fi [xiXj] ip [xgXt] + 



fc=l 
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